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Preface 


The  purpose  of  this  thesis  effort  was  to  begin  a  research  effort  into  the  develop¬ 
ment  of  Computer  Generated  Holography.  This  involved  the  development  of  a  system 
capable  of  producing  a  computer  generated  hologram. 

I  am  indebted  to  several  people  for  their  assistance  on  this  project.  First,  I  would 
like  to  thank  Major  Phil  Ambum,  my  thesis  advisor,  for  his  constant  encouragement, 
support,  and  advice.  Second,  I  would  like  to  thank  Dr.  Matt  Kabrisky  and  Major 
Steve  Rogers  who  served  on  my  thesis  committee.  Their  initial  efTorts  to  bring  me  up 
to  speed  on  holography  are  greatly  appreciated.  Next  I  would  like  to  thank  Captain 
Sean  Kelly,  AFSC  WRDC/KTD,  for  his  contribution  of  technical  information. 

Last,  I  wish  to  thank  my  wife  Debbie,  and  my  children  Karen  and  Benjamin  for 
their  constant  love  and  encouragement  during  the  last  18  months. 

Finally,  since  this  thesis  effort  deals  fundamentally  with  light,  and  since  I  wish 
to  give  credit  where  credit  is  due,  the  following  quote  is  appropriate. 

In  the  beginning,  God  created  the  heaven  and  the  earth.  And  the  earth 
was  without  form,  and  void;  and  darkness  was  upon  the  face  of  the  deep. 

And  the  spirit  of  God  move  upon  the  face  of  the  waters.  And  God  said.  Let 
there  be  light:  and  there  was  light.  And  God  saw  the  light,  that  it  was 
good:  And  God  divided  the  light  from  the  darkness. 

Genesis  1:1-4 
The  Authorized  Version 
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Abstract 

Computer  generated  holograms  are  used  for  a  variety  of  purposes.  One  promis¬ 
ing  approach  is  to  use  them  to  provide  three  dimensional  vic’vs  of  an  object.  The 
methods  currently  in  use  to  develop  computer  generated  holograms  use  the  Fast 
Fourier  Transform  and  are  not  geared  toward  developing  three  dimensional  images. 

This  thesis  effort  developed  a  method  to  produce  a  computer  generated  hologram 
by  implementing  an  equation  almost  identical  to  that  of  the  general  form  of  scalar 
diffraction  theory.  This  method  will  theoretically  allow  computer  generated  holograms 
to  be  made  for  a  wider  range  of  objects. 

The  interference  pattern  for  a  cube  was  calculated  using  this  method.  An 
electron-beam  lithography  machine  was  used  to  transform  this  pattern  to  a  glass 
plate. 


NON-FOURIER  COMPUTER 
GENERATED  HOLOGRAPHY  FOR  3-D 

DISPLAY 

I.  Introduction 

1. 1  Background 

Holography  is  a  technique  very  similar  to  photography.  The  primary  difference 
between  these  two  optical  recording  techniques  is  the  amount  of  information  that  is 
recorded.  Both  techniques  use  photographic  film  to  record  intensity.  Witn  photogra¬ 
phy  only  the  intensity  of  light  striking  the  photographic  film  is  recorded.  However, 
in  holography  the  intensity  pattern  recorded  on  the  film  has  phase  information  en¬ 
coded  in  it.  Because  holography  records  both  phase  and  intensity  information,  a 
hologram  can  be  used  to  produce  a  duplicate  of  the  wavefront  that  came  from  the 
object  recorded  in  the  hologram.  For  this  reason  holography  can  be  used  to  provide 
true  three-dimensional  imaging  while  photography  can  provide  only  two-dimensional 
imaging.  The  hologram’s  ability  to  produce  these  true  three-dimensional  images  has 
intrigued  people  for  many  years  and  has  prompted  much  research,  including  that 
reported  in  this  thesis. 

Holography  was  developed  in  1948  by  Gabor  in  an  attempt  to  reduce  the 
aberrations  of  electron  microscopes  (10).  Since  that  time  research  in  the  field  has 


equation  needed  to  produce  a  CGH  is 


H{X,,Yh)  =  J  J  J 


exp[i2ir\-\{X,  -  Xf^Y  +  -  YhY  +  (Z„  -  ZhY) 

((A'„  -  x,Y  +  {Yo  -  YhY  +  {Z,  -  z,Y)y^ 


2x1/21 


-dxodyodzo  (1) 


where  A  is  the  wavelength  of  light,  Xo,  Yo,  and  Zo  are  the  coordinates  of  a  point  on  the 
object,  and  Xh,  Yh  and  Zt,  are  the  coordinates  of  a  point  on  the  hologram  plane.  The 
integral  is  taken  over  the  entire  surface  of  the  object.  The  reader  unfamiliar  with  the 
mathematics  of  holography  is  referenced  to  Appendix  A. 

A  CGH  is  produced  by  using  equation  1  to  compute  the  intensity  value  at  a  finite 
number  of  points  in  the  hologram  plane.  These  intensity  values  are  then  rounded 
to  either  zero  or  one.  The  array  of  zeroes  and  ones  is  then  contoured  and  used  to 
drive  a  plotter.  The  resulting  interference  pattern  is  then  photographically  reduced 
to  the  appropriate  size  and  illuminated  with  a  smtable  light  source  to  reconstruct  the 
original  image. 

Computer  generated  holography  is  unique  because  the  object  of  which  a  hologram 
is  to  be  made  need  not  physically  exist.  Instead  only  a  mathematical  description 
of  the  object  is  needed.  As  a  result  computer  generated  holography  can  be  used  to 
produce  images  of  objects  that  would  be  very  difficult  or  impossible  to  produce  using 
standard  optical  techniques. 


1.2  Problem  Statement 

The  purpose  of  this  thesis  was  threefold.  First,  answer  the  following  question: 
Without  simplifying  equation  1,  is  the  calculation  of  a  CGH  a  tractable  problem  using 
the  computer  facilities  available  at  the  Air  Force  Institute  of  Technology  (AFIT)  and 
the  Air  Force  Weapons  Lab  (AFWL)  Super  Computer  Center?  Second,  develop  an 
initial  three-dimensional  CGH  capability  at  AFIT.  Finally,  begin  an  investigation 
to  determine  what  non-Fourier  simplifications  can  be  made  to  equation  1  without 
seriously  degrading  the  quality  of  the  resulting  holograms. 
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1.3  Current  Work 


Computer  generated  holography  ii  a  well  established  field  (20, 14,  22,  1).  Most  of 
the  work  in  this  area  deals  with  uses  of  holography  for  other  than  three-dimensional 
display.  Some  early  research  was  aimed  at  three-dimensional  display,  but  Hi  tie  of  the 
research  conducted  in  the  last  twenty  years  has  dealt  with  three-dimensional  display 
( 18).  Instead  recent  work  has  dealt  with  the  use  of  computer  generated  holography  in 
a  variety  of  other  fields  such  as  optical  data  processing,  interferometry,  information 
storage,  microscopy,  motion  pictures,  x-ray,  and  the  study  of  acoustics  (17,  23,  7,  13). 

1.4  Scope 

There  exist  many  methods  for  calculating  the  interference  pattern  of  a  CGH, 
Current  literature  contains  a  comprehensive  discussion  of  these  methods.  Chapter 
II  of  this  thesis  provides  an  overview  of  the  methods  found  in  the  literature.  The 
most  obvious  method  of  producing  a  CGH  is  a  two  step  process.  The  first  step  is  to 
use  equation  1  to  compute  the  values  at  a  finite  number  of  points  in  the  hologram 
plane.  The  second  step  involves  transferring  the  calculated  values  in  the  computer 
to  suitable  material  to  allow  illumination  and  hence  the  observation  of  the  original 
object.  The  computation  required  to  calculate  the  values  in  the  hologram  plane,  using 
this  method,  was  an  intractable  problem  when  computer  generated  holography  was 
first  considered. 

Because  the  obvious  and  straight  forward  way  of  producing  a  CGH  is  considered 
an  intractable  problem  research  has  been  directed  at  ways  of  simplifying  the  problem. 
One  simplification  which  seems  to  have  been  universally  implemented  is  to  reduce 
equation  1  into  a  Fourier  integral  and  then  use  the  Fast  Fourier  Transform  (FFT) 
algorithm  to  produce  results  (26:5).  It  appears  that  no  effort  has  been  made  to 
simplify  equation  1  without  converting  it  into  a  Fourier  integral.  Therefore,  this 
thesis  effort  concentrated  on  methods  of  simplification  which  did  not  involve  Fourier 
integrals  and  the  FFT. 
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Regardless  of  how  the  interference  pattern  is  calculated  the  pattern  must  be 
transferred  from  the  computer  onto  a  suitable  device  which  will  allow  illumination  of 
the  pattern.  There  are  three  primary  methods  available  to  transform  the  calculated 
values  into  an  actual  hologram.  Method  one  is  to  use  the  calculated  values  to  drive 
a  spatial  light  modulator  (16).  This  method  has  the  advantage  of  being  very  fast.  It 
has  the  disadvantage  of  being  temporary  because  when  the  spatial  light  modulator 
is  turned  off  the  hologram  is  gone.  Method  two  is  to  use  the  calculated  values  to 
drive  a  thermal  plastic  holographic  plate  (8).  This  method  is  slower  than  method  one 
and  has  the  same  disadvantage  of  being  only  temporary.  Method  three  is  to  draw 
the  calcxilated  interference  pattern  using  a  plotting  device  and  to  then,  if  necessary, 
photographically  reduce  the  plot  to  the  appropriate  size.  Method  three  is  by  far  the 
slowest  of  the  three  methods,  but  is  unique  because  it  produces  a  permanent  hologram 
that  can  be  reilluminated  and  viewed  multiple  times. 

The  primary  concern  of  this  thesis  effort  was  the  calculation  of  interference 
patterns  in  a  CGH  and  not  the  transformation  of  the  calculated  values  into  an 
actual  hologram.  For  this  reason  method  three  was  used  exclusively  to  transform 
the  calculated  values  into  the  actual  hologram.  This  allowed  multiple  holograms 
produced  by  different  methods  to  be  viewed  as  necessary  during  the  evaluation  phase 
of  the  research.  Figure  1  provides  an  illustration  of  the  research  conducted  for  this 
thesis. 

1.5  Approach 

This  thesis  effort  was  divided  into  four  distinct  steps:  a  literature  review, 
generation  of  interference  pattern  software,  generation  of  plotter  interface  software, 
and  an  evaluation  phase.  Each  of  these  steps  is  discussed  below. 

A  review  of  the  literature  on  computer  generated  holography  was  conducted. 
The  purpose  of  the  literature  review  was  two  fold.  First,  the  literature  review 
established  the  background  necessary  to  conduct  this  research  effort.  Second,  the 


literature  review  provided  an  understanding  of  the  current  methods  available  to 
produce  a  computer  generated  hologram.  Knowledge  of  current  methods  helped  to 
guide  the  direction  of  this  research. 

The  second  step  was  the  development  of  interference  pattern  software.  As 
mentioned  earlier  there  are  several  methods  available  to  accomplish  this  task.  Simple 
point  propagation  using  equation  1  was  chosen  as  the  first  method  to  implement. 
This  method  was  chosen  first  for  two  reasons.  First,  it  is  the  method  which  most 
accurately  simulates  the  development  of  a  classical  optical  hologram.  As  a  result, 
this  method  is  the  easiest  to  understand.  The  second  reason  for  choosing  this  method 
first  is  that  it  has  been  assumed  that  this  method  is  still  an  intractable  approach  to 
developing  CGH.  Given  the  advancement  in  computer  hardware  in  the  past  twenty 
years,  a  review  of  this  approach  was  warranted. 

The  third  step  was  to  develop  the  software  necessary  to  drive  the  e-beam 
lithography  machine.  This  step  involved  becoming  familiar  with  the  e-beam  input 
format  and  developing  an  algorithm  to  convert  the  calculated  intensity  values  into 
start  and  stop  points  for  the  e-beam  lithography  machine. 

The  final  step  in  the  process  was  the  evaluation  of  the  holograms  produced. 
This  evaluation  was  based  primarily  on  the  aesthetic  quality  of  the  images  produced. 
Image  brightness  and  clarity  were  the  major  factors  considered. 

1.6  Equipment  Required 

Successful  completion  of  this  thesis  required  the  following  equipment:  computers 
with  FORTRAN  and  C  compilers  for  the  calculation  of  the  interference  patterns,  a 
plotter  to  output  the  calc\ilated  interference  patterns,  and  a  gas  laser  to  illuminate 
the  hologram.  Computers  at  AFIT  and  AFWL  were  used  to  calculate  the  intensity 
values.  A  plotter  capable  of  writing  at  the  required  resolution  was  not  available  at 
AFIT,  so  the  e-beam  lithography  machine  at  the  Wright  Research  and  Development 
Center  (WRDC)  was  used.  Since  the  e-beam  lithography  machine  can  write  at  the 
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•  required  resolution,  no  photographic  reduction  was  required.  The  gas  lasers  used  to 

I  illuminate  the  holograms  were  available  at  AFIT. 

I 
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II.  Current  Knowledge 


2. 1  Introduction 

Since  1948,  when  Gabor  introduced  holography,  the  field  has  grown  to  include 
applications  in  a  variety  of  wavelengths  in  both  the  electomagnetic  and  acoustic 
ranges  (7,  13,  23).  Holograms  of  all  wavelengths  can  be  formed  digitally  and,  as  a 
result,  their  potential  was  recognized  early  on  and  investigation  into  their  use  began 
(18,  4,  24).  When  CGHs  are  made  using  optical  wavelengths  four  distinct  steps  are 
carried  out  (1:291): 

1.  Choosing  the  desired  object. 

2.  Calculating  the  object  wave  front  at  the  hologram  plane. 

3.  Encoding  this  wave  front  into  the  hologram  transmittance. 

4.  Optically  decoding  the  transmittance  by  illuminating  the  hologram  with  the 
reference  wave. 

CGHs  were  first  reported  by  Brown  and  Lohmann  in  1966  as  an  attempt 
to  develop  a  hologram  for  optical  spatial  filtering  (3).  Their  success  in  this  area 
stimulated  research  to  determine  applications  for  CGHs  and  to  determine  alternate 
methods  to  develop  CGHs.  An  important  aspect  of  computer  generated  holography 
was  also  noted  by  Brown  and  Lohmann.  The  problem  of  synthesizing  a  hologram 
is  the  opposite  from  typical  diffraction  problems.  With  normal  diffraction  problems 
the  object  is  given  and  the  diffracted  field  is  then  determined.  With  computer 
generated  holograms  the  image  is  prescribed  and  the  hologram,  or  diffracting  object, 
is  determined. 
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2. 2  Tripes  of  Computer  Generated  Holograms 

There  many  different  ways  to  digitally  produce  the  interference  patterns  neces¬ 
sary  to  make  a  hologram.  This  section  will  discuss  the  different  methods.  In  addition, 
the  work  currently  being  conducted  by  Steve  Benton  on  synthetic  holographic  stere¬ 
ograms  is  discussed. 

2.2.1  Actual  Computer  Generated  Holograms  The  methods  available  to  create 
true  computer  generated  holograms  can  be  divided  into  two  categories.  The  first 
category  is  point  propagation  and  involves  the  implementation  of  equation  1.  The 
second  category  involves  applying  two  simplifying  assumptions  to  equation  1  and 
uses  Fourier  Analysis  (26:4-5). 

2.2. 1.1  Equation  1  Implementation  Although  equation  1  is  the  fimda- 
mental  equation  used  to  mathematically  describe  the  intensity  patterns  produced  in 
a  hologram,  there  appears  to  be  no  published  accounts  of  its  direct  implementation. 
Presumably  this  is  due  to  the  fact  that  development  of  a  CGH  using  this  method  is 
considered  to  be  an  intractable  problem. 

2.2. 1.2  Fourier  Analysis  Methods  The  methods  discussed  in  this  section 
are  all  related  by  the  fact  that  Fourier  analysis  is  used  to  solve  the  problem.  The  use 
of  Fourier  analysis  is  desirable  because  this  allows  use  of  the  Fast  Fourier  Transform 
(FFT)  algorithm.  However,  before  the  problem  of  developing  a  CGH  can  be  properly 
synthesized  by  a  Fourier  Transform,  two  simplifying  assumptions  must  be  made. 
First,  if  the  dimensions  of  the  object  of  which  a  hologram  is  to  be  made  are  small 
in  comparison  with  the  distance  separating  the  object  and  the  hologram  plane,  then 
equation  1  can  be  reduced  to  a  Fresnel  integral.  Second,  if  the  hologram  plane  is 
small  in  comparison  with  the  distance  separating  the  object  and  the  hologram  plane, 
then  the  Fresnel  integral  resulting  from  the  first  assumption  can  be  further  reduced 
into  a  Fourier  integral  (26:4-5).  Thus,  if  these  two  assumptions  are  true,  the  use  of  the 
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FFT  is  a  valid  solution  to  the  integral  and  a  significant  time  savings  is  possible.  The 
remainder  of  the  section  is  devoted  to  discussion  of  different  approaches  to  developing 
CGHs  using  Fourier  analysis.  See  Appendix  B  for  more  detail  on  the  conversion  of 
equation  1  into  a  Fourier  integral. 

Cell  Oriented  Holograms  The  first  Computer  Generated  Holograms, 
developed  by  Brown  and  Lohmann,  were  originally  called  binary  holograms  (3).  This 
was  because  each  point  on  the  hologram  transparency  had  a  value  of  zero  or  one. 
That  is,  each  point  was  either  completely  transparent  or  completely  opaque.  Brown 
and  Lohmann  developed  their  CGHs  by  taking  the  Fourier  transform  of  an  image 
and  transferring  the  transform  values  to  the  hologram  transparency,  or  mask.  The 
transform  plane  was  divided  into  cells  of  equal  size.  Brown  and  Lohmann  discussed 
three  methods  of  representing  the  data.  In  the  first  method  each  cell  contained  two 
separate  apertures.  The  total  width  of  the  two  apertures  depended  on  the  magnitude 
of  the  transform  value  at  the  center  of  the  cell.  In  methods  two  and  three  each  cell 
had  only  one  aperture  whose  height  and  width  depended  on  the  magnitude  of  the 
transform  at  the  center  of  the  cell.  In  all  three  of  the  methods  the  lateral  position  of 
the  aperture  was  related  to  the  transform’s  phase  at  the  center  of  the  cell.  The  name 
detour  phase  hologram  was  later  applied  to  these  holograms  because  this  lateral 
shift  in  the  aperture  is  similar  to  diffraction  gratings  with  unequally  spaced  nilings 
(22:4353). 

Since  cell  oriented  CGHs  were  introduced  by  Brown  and  Lohmann  several 
different  variations  have  been  developed.  Dallas  reports  seven  distinct  types  of 
holograms  that  can  be  made  using  a  cell  oriented  approach  (1:295-311).  They  are 
all  similar  and  are  based  on  the  detour  phase  holograms  developed  by  Brown  and 
Lohmann.  Tricoles  provides  the  following  excellent  explanation  of  the  equations  used 
for  detour  phase  holograms.  Figure  2  explains  several  of  the  variables  used  in  the 
following  equations. 


The  aim  is  a  hologram  with  transmittance  H{vx,Vy),  which  produces  an 
image  amplitude  u(x,y)  when  illuminated  by  a  plane  wave  exp(z27rj/3.io)  = 
E{v^xq)  in  the  arrangement  of  figure  3,  where  xh  =  A/i/,  and  yh  =  A/i/„. 
The  requirement  is  that  h{x,y)  the  amplitude  diffracted  by  the  hologram 
be  proportional  to  that  of  the  image  in  a  finite  region;  that  is, 

h{x,y)  =  Ku{x,y).  (2) 


I'x 


Figure  2.  TVpical  cell  in  a  binary  hologram.  Position  indices  are  n  and  m.  W^dth  is 
CSu,  height  VF„„5i/,and  lateral  displacement  from  the  cell  center  is  PnmSv 
(22) 


With  the  image  width  Ax  and  height  Ay,  the  diffracted  amplitude  is 


11 


Figure  3.  Optical  setup  for  reconstruction  with  a  point  source  (SOU)  at  lo-  HOLO, 
IMG,  and  FREQ  are,  respectively,  the  hologram,  image  and  frequency 
planes. (22:4353) 


h{x,y)  =  rect{x/Ax)rect{y/Ay)J  J  H{i/x,Uy)E[{x  +  xo)ug  +  yi/y]duxduy,  (3) 

where  rectz  has  a  unit  value  for  \z\  <  1/2  and  zero  otherwise. 

The  hologram  is  represented  by  sampled  values.  The  number  of  samples 
depends  on  image  size  and  image  resolution.  If  6x  is  the  resolution  element 
size  in  the  image,  the  number  of  resolvable  points  in  the  image  is,  with 
Sz  = 

=  {A,A,)nS.f  =  (A.A,)(A:^f.  (4) 

N  is  called  the  number  of  degrees  of  freedom.  The  hologram  is  assumed  to 
have  at  least  points  to  preserve  information;  thus 

>  (Ai/5x)^  =  (AxAt/)^.  (5) 

In  fact,  Su  was  taken  to  be  (Ax)"^  on  the  basis  of  the  sampling  theorem. 

To  develop  expressions  for  aperture  dimensions,  consider  equation  (  2) 
and  apply  scalar  diffraction  theory  to  the  left-hand  side.  The  hologram’s 
transmittance  is 

H{v^Uy)  =  -  (n  +  Pnm)Su]lc8v  X  rect[{vy  -  m6v)IWj,rnSv].  (6) 

n  m 

The  rectangle  functions  describe  the  aperture  dimensions  in  a  cell  as  in 
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(7) 


figure  3.  Illumination  by  a  plane  wave  produces  complex  amplitude 

^  =  J  J  H{l/j:Uy)E[{x  +  xo)u^  +  yi/y]dl/j;dUy 

=  c{6u)^3inc[cSu{x  +  xo)]  X]  Y1  ^ 

sinc{yWr,mSi/)E{Si/[{x  +  xo)(n  +  P„m)  +  yrn]}  (8) 

For  the  right-hand  side  of  equation  (  2)  the  described  image  u(x,y)  in  its 
Fourier  representation  with  the  tilde  signifying  the  transform  is 

u(x,y)  =  J  j  u  {v^Uy)E{xVj:  -I-  yUy)dv^dvy.  (9) 

Consider  the  Fourier  sampling  theorem 

u=  ^  X]  ^  {nlAx,mf  Ax)  x  sinc{i/xAx  —  n)sinc{vyAx  —  m)  (10) 

n  m 

With  equation  (  10),  equation  (  9)  gives 

u(x,y)  =  rect(xlAx)  x  rect{y  j  Ax)  ^  X^  u(n^i/,  m(5)F[(5i/(xn -f  ym)].  (11) 

n  m 

To  satisfy  the  condition  in  equation  (  2),  equation  (  8)  is  approximated  to 
equal  equation  (  11)  term  by  term.  The  result  is 

{6u)'^WnmE[xQ8u{n  +  P„m)]  ^  const,  u  {n6v,m6i/),  (12) 

where  const,  u  {n6i/,m6u)  is  defined  as  c(^i/)^A„„F(0„„/27r),  and  and 
4)nm  are  the  magnitude  and  phase  of  the  transform  u.  Thus 

(13) 

Pnm  -I-  n  «  <l>„m/2Trxo6i/.  (  14) 

In  words,  controls  magnitude;P„m  controls  phase.  If  xq  is  chosen  so 
that  xqSv  is  an  integer  M, 


Pnm^<i>nm/27rM  (15) 

Equation  (  15)  shows  that  the  lateral  displacement  is  proportional  to  the 
phase  of  the  transform,  an  important  relation,  which  has  intuitive  meaning 
as  the  phase  difference  of  waves  originating  from  two  separated  Huygens 
sources  (22:4353-4354). 


Point  Oriented  Holograms  For  point  oriented  holograms  it  is  re¬ 
quired  that  the  complex  amplitude  to  be  recorded  in  the  hologram  be  known  for  each 
point  on  the  object  (1:311).  The  complex  amplitude  of  each  point  on  the  object  is  then 
mapped  to  each  point  on  the  hologram  plane.  This  method  has  the  potential  of  being 
very  accurate  because  it  uses  as  its  basis  ^’he  equations  of  point  light  propagation 
as  discussed  in  Appendix  A.  An  exact  approximation  of  an  optical  hologram  can  be 
theoretically  achieved  using  this  method  if  an  infinite  number  of  points  are  used 
to  describe  an  object  and  an  infinite  number  of  points  are  allowed  in  the  hologram 
plane.  This  of  course  is  not  practical  so  several  variations  to  this  method  have  been 
developed  (1,  18:311-320). 

Plane  Propagation  The  methods  discussed  thus  far  have  been  ap¬ 
proximations  to  the  ’exact’  method  of  calculating  the  transmittance  values  present  at 
the  hologram  plane.  These  approximations  have  been  developed  because  the  ’exact’ 
method,  which  specifies  all  depth  information,  is  considered  computationally  expen¬ 
sive.  This  section  is  devoted  to  a  discussion  of  the  ’exact’  method  and  an  algorithm 
used  to  implement  this  method  (1:337-338). 

First,  consider  the  reconstruction  image  shown  in  fig.  4.  The  image  found  in  the 
image  plane,  zq,  is  planar  and  therefore  two-dimensional.  It  is  possible  to  move  the 
image  a  distance  Az  out  of  the  image  plane  by  building  into  the  CGH  an  auxiliary  lens 
of  the  proper  focal  length.  Such  an  auxiliary  lens  can  be  described  mathematically  as 

/e  =  [1/(20  -I-  Az)]"^  =  [zo(zo  +  Az)]/(5z 

This  lens  can  actually  be  incorporated  into  the  CGH  by  multiplying  the  object 
Fourier  spectrum  by  the  quadratic  phase  factor 

eip[t7r(x^-f-y^)/A/e] 
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Figure  4.  Fourier  hologram  reconstruction  setup  (?:291) 

Rewriting  the  quadratic  phase  factor  in  terms  of  the  reduced  coordinate  system  yields 

exp{iTrXAz[{zo  +  Az)/zo\{fi^  -h  i/^)} 

The  Fourier  spectrum  for  the  displaced  object  becomes 

uj  (/i,  I/)  =u  (p,  i/)exp[t7rXAz(/x^  + 

This  equation  represents  the  translation  of  a  plane  by  Az  along  the  optical  axis. 

If  an  object  is  described  as  a  collection  of  planes  this  translation  can  be  used  to 
propagate  each  plane  from  its  original  position  to  the  hologram  plane.  As  each  plane 
is  translated  the  complex  amplitude  is  summed  and  multiplied  by  the  transmittance 
of  the  plane  currently  being  translated.  This  method  is  explained  in  detail  by  Ichioka, 
Izumi,  and  Suzuki  (15).  Dallas  provides  the  following  analysis  of  this  method. 

Because  the  multiplication  of  the  incident  wave  by  the  plane’s  transmit- 
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tance  is  performed  in  the  object  domain  and  the  propagating  takes  place 
in  the  Fourier  domain,  one  is  constantly  transforming  from  one  domain  to 
the  other;  hence  the  term  ping-pong  propagation.  K  there  are  Np^  planes 
in  the  object  then  2Npl  -  1  DFTs  are  performed  and  Npl  ~  1  transmittance 
multiplications.  One  must  be  very  careful  that  the  space-bandwidth  prod¬ 
uct  of  the  hologram,  number  of  resolution  cells,  is  large  enough.  The  image 
space  in  the  last  plane  must  be  surroimded  by  a  zero  frame  which  is  large 
enough  to  essentially  capture  waves  expanding  from  previous  planes.  A 
good  rule  of  thumb  is  to  multiply  the  reciprocal  of  the  CGH’s  /-number  at 
the  first  plane  by  the  distance  from  the  plane  of  interest  to  the  final  plane. 
Dividing  by  the  sampling  distance  in  the  final  plane  gives  the  number  of 
samples  in  the  zero  frame.  (1:338) 

2.2.2  Synthetic  Holographic  Stereograms  Benton  has  conducted  the  only  re¬ 
cent  research  aimed  at  using  CGHs  for  three-dimensional  display.  Benton’s  approach 
is  as  follows: 

A  series  of  perspective  views,  corresponding  to  the  expected  viewpoints  of 
the  audience,  are  computed  by  conventional  computer-graphic  techniques. 

Each  perspective  view  is  then  projected  with  laser  light  onto  a  sheet  of  high 
resolution  film  from  the  angle  corresponding  to  its  computed  viewpoint, 
while  it  is  overlapped  by  a  second  coherent  "reference”  beam  to  produce  a 
holographic  exposure  that  records  the  direction  of  the  image  light.  After 
all  the  views  have  been  recorded  in  this  way,  the  hologram  is  processed 
and  then  illuminated  so  that  each  view  is  sent  back  out  in  the  direction  it 
was  projected  from.  That  is,  toward  its  intended  viewing  location,  so  that 
a  viewer  moving  from  side  to  side  sees  a  progression  of  views  as  though 
he  of  she  were  moving  around  an  actual  object.  If  the  perspectives  are 
accurately  computed  and  registered,  the  resulting  composite  image  will 
look  like  a  solid  3-D  object.  Such  a  composite  or  synthetic  hologram  is 
termed  a  "holographic  stereogram.”  It  mimics  the  visual  properties  of  a  true 
hologram,  but  lacks  its  enormous  information  content  and  interferometric 
accuracy.  (2) 


2.3  Conclusion 

The  potential  benefits  of  computer  generated  holography  are  well  known  and  as 
a  result  much  research  has  gone  into  developing  viable  methods  of  producing  CGHs. 
Most  of  the  work  in  this  area  has  been  devoted  to  developing  methods  applicable 
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to  disciplines  other  than  three-dimensional  display.  There  has,  however,  been  some 
success  in  developing  CGHs  for  three-dimensional  display  (18,  1).  It  appears  that  all 
work  done  so  far  in  the  area  of  three-dimensional  CGH  has  used  some  form  of  Fourier 
analysis.  The  use  of  Fourier  analysis  places  a  restriction  on  the  size  of  the  object 
and  the  size  of  the  resulting  hologram.  It  is  possible  that  the  tremendous  growth  in 
computing  power  over  the  past  20  years  wall  now  allow  research  to  begin  on  methods 
of  developing  CGHs  without  the  uses  of  Fourier  analysis. 
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III.  System  Development 


3. 1  Introduction 

This  chapter  details  the  development  of  the  software  system  built  during  this 
thesis  effort.  In  addition,  the  different  hardware  configurations  used  throughout  the 
thesis  effort  are  also  discussed  in  this  chapter. 

The  software  development  was  divided  into  two  categories.  The  first  was 
development  of  the  software  to  calculate  the  intensity  values  present  in  the  hologram 
plane.  The  second  was  development  of  the  software  to  generate  the  Calma  stream 
format  to  drive  the  e-beam  lithography  machine. 

All  software  was  originally  written  in  the  C  programming  language  and  was 
run  on  a  Sun  4  and  a  DEC  GPX.  Subsequent  versions  of  all  software  were  run  on  the 
Cray  2  located  at  the  Supercomputer  Center,  Kirtland  AFB,  NM.  The  unreliability 
of  the  Cray  C  compiler  resulted  in  all  software  being  converted  to  FORTRAN  before 
reliable  results  could  be  obtained  from  the  Cray. 

3.2  Intensity  Value  Software 

The  Intensity  Value  Software  (IVS)  went  through  three  distinct  phases  of 
development.  Phase  one  produced  Run  Length  Encoded  (RLE)  files  which  allowed  the 
resulting  interference  patterns  to  be  observed  on  the  screen  of  a  graphics  workstation. 
Phase  two  produced  output  in  the  form  of  a  file  containing  a  value  between  0  and  255 
for  each  of  the  points  in  the  hologram  plane.  Phase  three  consisted  of  translating  the 
phase  two  program  written  in  C  into  FORTRAN  to  allow  maximum  optimization  and 
vectorization  on  the  Cray  2. 

All  three  phases  of  the  IVS  had  the  same  general  function.  That  function  was 
to  implement  equation  1  as  modified  in  Appendix  A,  section  A.4.  This  was  done  by 
keeping  a  real  and  an  imaginary  part  for  each  element  in  the  array  that  served  as 
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the  hologram  plane.  For  each  point  on  the  hologram  plane,  the  real  and  imaginary 
contribution  from  each  point  on  the  object  was  calculated  using  equation  25.  After  the 
contribution  from  each  point  had  been  calculated,  the  contribution  from  the  reference 
beam  was  calculated  for  every  point  on  the  hologram  plane.  At  this  point  the  complex 
conjugate  of  the  resisting  complex  numbers  was  taken  in  order  to  produce  the  overall 
intensity  at  each  point  on  the  hologram  plane. 

Equation  25  was  approximated  by  using  three  loops.  The  outer  most  loop  was 
executed  once  for  every  point  on  the  object.  The  second  loop  was  executed  once  for 
every  X  value  in  the  array  used  as  the  hologram  plane.  The  inner  loop  was  executed 
once  for  every  Y  value  in  the  array  used  as  the  hologram  plane.  The  only  difference 
between  phase  one  and  phase  two  IVS  was  the  form  in  which  the  resulting  intensity 
values  were  saved.  Phase  three  was  identical  to  phase  two  except  that  phase  three 
was  coded  in  FORTRAN,  not  C. 

3.2.1  Phase  One:  RLE  Output  The  purpose  of  phase  one  IVS  was  to  insure 
equation  25  had  been  properly  programmed.  This  was  done  by  calculating  the 
interference  patterns  for  a  few  very  simple  objects  for  which  the  interference  patterns 
were  easily  predictable.  The  hologram  plane  consisted  of  a  300  X  300  array  of  values. 
The  calculated  values  were  normalized  by  dividing  each  of  the  values  by  the  largest 
value.  Since  the  Sun  4  is  capable  of  providing  256  shades  of  gray,  each  of  the 
normalized  values  was  scaled  to  between  0  and  255.  A  linear  scale  was  used.  The 
300  X  300  array  of  scaled  values  was  then  converted  into  a  RLE  file  for  viewing. 

The  first  images  to  be  rendered  using  phase  one  software  were  the  interference 
patterns  created  by  two  point  light  radiators  located  at  distances  of  10,  25,  50,  and 
100  wavelengths  from  the  hologram  plane.  A  photograph  of  a  composite  of  these 
four  images  is  shown  in  figure  5.  Clockwise  from  the  upper  left,  the  image  shows 
the  interference  pattern  formed  when  the  separation  between  the  points  and  the 
hologram  plane  is  10,  25,  50,  and  100  wavelengths  respectively. 
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Based  on  these  images,  it  was  determined  that  equation  25  had  been  properly 
implemented  and  work  began  on  phase  two  IVS. 

3.2.2  Phase  Two:  Binary  Output  Phase  two  IVS  was  identical  to  phase  one 
IVS  in  its  calculation  of  the  hologram  plane  intensity  values.  Only  the  output  method 
was  different.  With  phase  one  software,  the  scaled  values  between  0  and  255  were 
used  to  produce  an  RLE  image.  With  phase  two  software,  the  scaled  values  were 
simply  output  to  a  file.  These  scaled  values  were  later  used  by  the  plotter  software  to 
drive  the  e-beam  lithography  machine. 

Phase  two  IVS  was  uploaded  to  the  Cray  2  to  determine  what  speed  improve¬ 
ments  could  be  realized  by  running  on  the  Cray.  It  was  discovered  that  the  Cray  2 
C  compiler  was  unreliable.  The  major  problem  found  was  that  the  SIN  and  COS 
functions  do  not  return  correct  values  for  large  numbers.  Since  the  use  of  the  SIN  and 
COS  functions  were  vital  to  the  correct  operation  of  the  program,  it  was  determined 
that  the  final  version  of  the  program  should  be  coded  in  FORTRAN. 

3.2.3  Phase  Three:  Conversion  to  FORTRAN  Phase  three  IVS  was  a  line  by 
line  conversion  of  phase  two  IVS  into  FORTRAN.  In  phase  three  IVS  the  hologram 
plane  was  represented  by  a  36,000  X  36,000  element  array.  Each  element  of  the 
array  represented  a  point  on  the  actual  hologram  plane  that  was  separated  by  .5 
microns  from  each  of  its  adjacent  elements.  This  means  that  the  array  represented 
a  hologram  measuring  1.8cm  X  1.8cm.  Since  no  computer  available  had  sufficient 
memory  to  hold  a  36,000  X  36,000  element  array,  the  values  in  the  hologram  plane 
were  iteratively  caloilated  in  2000  X  2000  arrays.  This  resulted  in  the  large  array 
being  subdivided  into  324  smaller  arrays,  each  of  which  was  labeled  by  its  row  and 
column  letter.  Figure  6  shows  the  layout  of  the  large  array  and  the  labeling  of  the 
subarrays. 
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Figure  6.  Hologram  plane  array  and  its  subdivisions 


3.3  Plotter  Software 


3.3.1  Calma  Electron-Beam  Machine  Once  the  IVS  was  completed  and  had 
been  executed,  it  was  necessary  to  write  a  program  to  interface  the  e-beam  lithography 
machine  to  the  data  produced  by  the  IVS.  The  e-beam  machine  used  was  produced 
by  General  Electric  under  the  name  Calma.  The  e-beam  machine  is  part  of  a  larger 
computer  aided  design  (CAD)  system  used  to  develop  very  large  scale  integration 
(VLSI)  circuitry.  A  computer  workstation  is  used  to  design  the  required  circuitry. 
The  designed  circuit  is  then  converted  into  the  Calma  GDSIF’^  stream  format  by 
the  workstation.  The  tape  containing  the  circuit  design,  in  stream  format,  is  then 
provided  to  the  e-beam  lithography  machine  which  etches  the  designed  circuit  into  a 
quartz  plate.  One  of  the  advantages  of  this  system  is  that  a  tape  containing  data  in 
the  Calma  GDSII^^  stream  format  can  be  read  back  into  the  workstation  and  viewed 
on  the  screen. 

The  hologram  produced  by  the  e-beam  machine  was  actually  a  quartz  plate 
coated  with  chrome.  The  e-beam  machine  produced  the  required  interference  pattern 
by  removing  the  chrome  at  only  the  locations  of  high  intensity.  The  areas  where  the 
chrome  has  been  removed  are  interference  maxima  and  the  areas  where  the  chrome 
has  not  been  removed  are  interference  minima.  See  Appendix  A,  section  A.2  for  an 
explanation  of  interference  m2udma  and  minima.  The  use  of  an  e-beam  lithography 
machine  was  of  significant  benefit  because  no  photographic  reduction  was  required. 
This  was  due  to  the  submicron  writing  resolution  capability  of  the  e-beam  machine. 

3.3.2  Intensity  Value  to  Calma  Stream  Format  Interface  Software  The  first 
plotter  software  written  was  an  attempt  to  become  familiar  with  the  Calma  GDSII^^ 
stream  format.  The  program  simply  produced  the  Calma  commands  to  draw  two 
perpendicular  lines.  The  tape  containing  these  commands  was  taken  to  WRDC 
and  loaded  into  the  Calma  workstation  for  viewing.  Subsequent  version  of  this 
first  program  produced  the  Calma  commands  to  draw  increasingly  complex  patterns, 


including  lines  of  varying  width.  Each  version  was  viewed  on  the  Calma  workstation 
to  insure  the  desired  results  had  been  attained.  This  exercise  provided  the  necessary 
background  to  develop  the  actual  interface  to  the  e-beam  lithography  machine. 

The  software  which  interfaced  the  calculated  intensity  values  to  the  e-beam 
lithography  machine  functioned  in  the  following  way.  The  user  supplied  the  name 
and  location  of  one  of  the  small  arrays  in  figme  6  to  the  program.  As  each  array  was 
supplied  to  the  program,  a  threshold  was  used  to  convert  the  intensity  values  into  a 
binary  form.  All  values  above  the  threshold  were  considered  on  and  all  values  below 
the  threshold  were  considered  off.  Next,  the  program  began  at  the  bottom  of  the  first 
column  and  searched  upward  imtil  it  foimd  the  first  value  above  the  threshold.  The 
location  of  this  value  was  stored  and  the  program  continued  searching  upward  imtil 
a  value  below  the  threshold  was  fovmd.  The  location  previous  to  the  one  below  the 
threshold  was  stored  and  the  program  produced  the  stream  commands  necessary  to 
produce  a  line  between  the  two  stored  locations.  If  the  beginning  and  ending  points 
were  the  same,  the  program  looked  to  the  right  and  left  to  determine  if  this  point  was 
part  of  a  horizontal  line  of  length  greater  than  one.  If  it  was,  the  point  was  not  drawn 
at  that  time.  The  program  searched  every  column  in  this  manner  until  the  entire 
array  was  complete.  At  that  point  the  program  began  with  the  first  row  searching 
for  horizontal  lines  in  the  same  manner.  Any  vertical  line  of  length  one  which  was 
skipped  in  the  earlier  pass  was  formed  at  this  time.  This  two  pass  approach  was 
used  to  reduce  the  total  size  of  the  stream  format  files.  If  the  interference  pattern 
has  many  horizontal  lines  in  it,  a  one  pass  method  using  a  vertical  search  would  be 
required  to  write  the  stream  commands  for  many  paths  of  length  one.  However,  by 
skipping  paths  of  length  one  when  they  are  part  of  a  horizontal  path  reduces  the  total 
number  of  commands  which  must  be  produced. 

The  original  plan  was  to  transport  the  data  resulting  from  the  IVS  from  Kirtland 
AFB,  NM  to  AFIT.  However,  transporting  1.3  gigabytes  of  data  turned  out  to  be  a 
difficult  task.  As  a  result,  the  plotter  software  which  was  originally  written  in  C, 
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was  transported  to  the  Cray  2.  Again  it  was  discovered  the  Cray  C  compiler  was 
unreliable.  For  instance,  while  trying  to  write  out  the  Calma  commands  to  a  disk  file 
it  was  discovered  that  the  decimal  number  10  could  not  be  written  out.  Instead  some 
other  value  was  written  to  the  file.  So  this  program  was  also  converted  from  C  into 
FORTRAN. 
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TV.  Results  I  Conclusion 


4. 1  Introduction 

This  chapter  contains  a  discussion  of  the  results  and  conclusions  drawn  from 
this  thesis  effort.  There  is  a  section  on  the  intensity  value  software,  the  plotter 
software,  the  hologram  produced,  and  the  conclusions  drawn. 

4.2  Results 

4.2.1  Intensity  Value  Software  The  IVS  was  used  to  calculate  the  object  wave 
front  present  at  a  hologram  plane.  The  object,  in  this  case,  was  a  cube.  Figure  7 
shows  the  relationship  of  the  cube  to  the  hologram  plane.  The  hologram  plane  and 
face  A  of  the  cube  are  parallel  to  one  another  and  are  separated  by  50,000  microns. 
The  cube  measures  50,000  X  50,000  X  50,000  microns.  Beginning  in  the  upper 
left,  moving  clockwise,  the  coordinates  of  the  comers  of  face  A  are  (-50000,50000,0), 
(0,50000,0),  (0,0,0),  and  (-50000,0,0).  The  top,  back  comer  of  face  B  is  at  coordinate 
(0,50000,50000)  and  the  bottom,  back  comer  of  face  B  is  at  coordinate  (0,0,50000).  All 
coordinates  are  given  in  m'  rons.  A  total  of  55  points  were  used  to  define  the  edges  of 
the  cube.  Because  of  t^  orientation  of  the  hologram  plane  and  its  size,  only  faces  A 
and  B  of  the  cube  are  visible  at  me  hologram  plane.  As  a  result,  all  55  points  defining 
the  cube  lie  on  the  7  edges  defining  these  2  faces. 

A  36,000  X  36,000  element  array  represented  the  hologram  plane.  Each  element 
of  the  array  contained  a  real  and  an  imaginary  part.  During  the  calcrJation,  the 
real  and  imaginary  contribution  from  each  point  on  the  object  was  added  to  the 
appropriate  part  of  each  element  of  the  array.  After  the  intensity  values  in  the  first 
section  of  the  hologram  plane  were  calculated,  the  IVS  was  stopped.  A  program  was 
written  to  convert  the  calculated  intensity  values  into  and  RLE  image.  This  was  done 
in  an  effort  to  detect  any  errors  that  may  have  occiured  in  the  calculation.  Figure  8 
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is  a  photograph  of  the  resulting  RLE  image.  This  image  was  determined  to  be  correct 
and  the  IVS  was  restarted  and  allowed  to  calculate  the  intensity  values  for  the  entire 
array.  After  the  calculation  was  complete,  the  complex  conjugate  was  taken  at  each 
element  of  the  array.  The  resulting  real  number  was  then  stored  as  the  intensity  at 
that  point  on  the  hologram  plane.  After  the  normalization  and  scaling,  as  discussed 
in  chapter  III,  section  3.2.1,  the  resvdting  one  byte  for  each  of  the  1,296,000,000 
elements  of  the  array  was  written  out  to  a  disk  file.  Calculation  of  this  1.3  gigabytes 
of  data  took  approximately  59.5  hours  of  Cray  2  CPU  time. 


Figure  8.  Photograph  of  Cube  Interference  Pattern 

4.2.2  Plotter  Software  Each  of  the  small  arrays  shown  in  figure  6,  which 
contained  the  intensity  values  for  that  part  of  the  hologram  plane,  was  supplied  to 
the  plotter  software  which  created  the  appropriate  Calma  commands  to  produce  the 
interference  fringes  contained  in  that  array.  The  calciilated  intensity  values  and  the 


28 


required  interference  fringes  were  related  as  follows.  Since  the  e-beam  machine  had 
no  gray  scale  capability,  a  threshold  was  used  to  convert  numbers  into  binary  form. 
All  values  above  the  threshold  were  considered  on,  or  of  high  intensity  and  all  values 
below  the  threshold  were  considered  off,  or  of  low  intensity.  Therefore,  if  several 
adjacent  elements  of  the  array  were  above  the  threshold,  an  interference  maxima 
was  formed  there.  Likewise,  if  several  adjacent  elements  of  the  array  were  below  the 
threshold,  an  interference  minima  was  formed  there. 

Production  of  all  the  commands  necessary  to  draw  the  hologram  on  the  e-beam 
machine  required  approximately  5  hours  of  Cray  2  CPU  time. 

4.2.3  The  Hologram  The  Calma  commands  produced  by  the  plotter  software 
did  not  directly  interface  to  the  e-beam  machine.  Instead,  the  commands  for  each 
section  of  figure  6  had  to  be  converted  to  a  format  which  was  directly  readable  by 
the  e-beam  machine.  No  documentation  was  available  for  this  second  format,  so  the 
plotter  software  could  not  write  this  format  directly.  At  the  time  of  this  writing  the 
conversion  process  was  approximately  50%  complete.  Each  file  containing  commands 
for  one  section  of  the  hologram  plane  takes  approximately  45  minutes  to  convert. 
The  e-beam  machine  operators  estimate  the  e-beam  machine  will  require  48  hours  to 
draw  the  entire  pattern  once  the  conversion  is  complete. 

4. 3  Concl  usions 

This  thesis  effort  had  three  primary  purposes  as  stated  in  Chapter  1,  section  1.2. 
Each  of  these  purposes  is  discussed  below. 

The  first  and  most  important  purpose  of  this  thesis  was  to  determine  if  the 
calculation  of  the  information  necessary  to  make  a  computer  generated  hologram  of 
an  object  was  a  tractable  problem  using  equation  1.  To  answer  this  question  the 
calculation  was  carried  out  for  a  very  simple  geometric  object.  Computers  in  the  class 
that  are  typically  available,  such  as  a  Sun  4  or  a  Vax  11/780,  proved  to  be  incapable 
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of  making  this  calculation.  This  was  due  to  the  extremely  large  numbers  used  during 
the  calcvilation.  Machines  in  the  class  of  a  Cray  were  required  in  order  to  obtain 
accurate  results.  A  Cray  2  was  used  to  make  the  calculation  and  required  almost 
60  hours  of  CPU  time  to  complete  the  calculation.  The  result  of  this  calculation 
was  1.3  gigabytes  of  data.  Most  computer  facilities  do  not  have  1.3  gigabytes  of  disk 
space  available  to  a  single  user.  In  light  of  the  excessive  time  required  and  the  large 
amount  of  data  produced,  it  is  my  opinion  that  this  calculation,  as  implemented,  is 
not  a  tractable  problem.  This  method  can  be  used  only  at  facilities  with  access  to  the 
equipment  mentioned  above. 

The  second  purpose  of  the  thesis  was  to  develop  an  initial  three  dimensional 
computer  generated  holography  capability  at  AFIT.  Fiilfillment  of  this  goal  required 
the  ability  to  make  the  necessary  calculations  and  to  convert  the  calculated  intensity 
values  into  a  form  suitable  for  viewing  the  three  dimensional  object.  This  goal  was 
met.  As  stated  above,  a  Cray  2  was  used  to  calculate  the  intensity  values.  These 
values  were  then  used  to  drive  a  program  which  produced  commands  for  an  e-beam 
lithography  machine.  The  e-beam  lithography  machine  then  drew  the  interference 
pattern  on  a  chrome  covered  quartz  plate. 

The  final  purpose  of  this  thesis  effort  was  to  begin  an  investigation  into  non- 
Fourier  simplifications  that  could  be  made  to  equation  1.  This  investigation  revealed 
two  approaches  which  might  be  pursued. 

1.  The  first  avenue  which  could  be  investigated  is  an  obvious  one  which  is  often 
applied  to  difficult  integral  equations.  This  method  is  to  use  Monte  Carlo 
techniques  to  approxim.ate  the  solution  to  the  equation.  There  are  multiple 
forms  of  these  techniques.  If  one  chooses  to  investigate  this  avenue  further 
care  must  be  taken  to  insure  the  proper  techniques  are  applied  to  this  equation. 
Several  books  have  been  written  about  Monte  Carlo  techniques  and  their 
applicability  to  different  types  of  equations.  One  particularly  good  book  is  The 
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Monte  Carlo  Method  by  Sobol  (21). 


2.  The  second  possible  area  for  further  research  deals  with  making  estimations  for 
part  of  the  equation.  As  discussed  in  appendix  A,  section  A.4  equation  1  can  be 
altered  to 


cosd  +  isinO 

{{Xo  -  Xh)^  +  (Vc  -  n)2  +  (Z,  -  Z;.)2)l/2 


where  ff  =  2rr)<-\(X,  ~  X^f  +  +  (Z,  - 

Preliminary  investigation  suggest  that  an  accurate  estimate  for  the  value  of 
6  may  be  obtained  by  using  a  step  function.  In  calculating  6,  part  of  the 
equation  is  constant,  namely  2TrX~^.  The  remainder  of  the  equation  is  the 
distance  separating  a  particular  point  on  the  object  from  a  particular  point  on 
the  hologram  plane.  This  distance  is  calculated  by  substituting  the  {X,  Y,  Z) 
coordinate  of  a  point  on  the  object  for  the  variables  Xo,Yo,  and  Zo  respectively. 
Then  the  (X,  Y,  Z)  coordinate  of  a  point  on  the  hologram  plane  is  substituted  for 
the  variables  X/,,  Yh,  and  Zh  respectively.  The  appropriate  arithmetic  is  carried 
out  and  the  result  is  multiplied  by  the  constant  mentioned  above.  The  result 
is  the  value  of  6.  The  value  of  9  could  be  estimated  if  the  distance  part  of  this 
calculation  could  be  accurately  estimated.  This  would  result  in  a  sigmficant 
time  savings. 

Figure  9  shows  the  geometry  of  this  problem  and  is  a  side  view.  Figure  10  is  a 
front  view  of  the  same  geometry.  As  can  be  seen  from  figure  9,  the  distance  being 
calculated  is  the  length  of  the  hypotenuse  of  a  right  triangle.  In  figure  9,  when 
9  is  calculated  for  each  point,  Pi  -  P4,  all  values  in  the  equation  are  the  same 
except  the  Y  value  in  the  hologram  plane.  For  each  successive  calculation  of  9, 
this  value  is  simply  incremented  by  one.  If  the  values  (X^  -  X/,)^,  {Yo  -  Yhf,  and 
(Zo  -  Z^y  are  large,  incrementing  Yh  by  one  will  have  very  little  effect.  Therefore, 
given  the  distance  from  a  point  on  the  object  to  a  point  on  the  hologram  plane. 
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the  distance  to  the  next  point  on  the  hologram  plane  should  be  very  nearly  the 
same  and  predictable.  If  this  is  the  case,  successive  values  of  9  should  also  be 
very  predictable.  If  this  hypothesis  holds  true  and  the  value  of  each  successive 
9  is  predictable  and  changes  only  a  very  small  amount,  it  may  be  possible  to 
estimate  the  am  and  cos  of  each  successive  9. 


Hologram  Plane 


{Xh,Yh  +  2,Zh) 
{Xh,Yh  +  1,Zs) 
{Xh,Yk,Zk) 


Figure  9.  Side  View  of  the  Problem  Geometry 


Figure  10.  Front  View  of  the  Problem  Geometry 
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Appendix  A.  Holography 


A.  1  Introduction 

This  appendix  contains  three  parts.  Part  one  is  a  general  discussion  of  the  topic 
of  holography  and  other  relevant  topics  in  optics.  It  contains  historical  informa*  ”> 
as  well  as  a  description  of  the  general  methods  used  to  produce  an  optical  hologram. 
Part  two  contains  a  presentation  of  the  mathematical  theory  behind  holography.  Part 
three  contains  a  discussion  of  the  calculation  of  Uo- 

Many  authors  have  written  books  directed  specifically  at  holography  and  its 
uses  (5,  6,  12).  In  addition,  most  contemporary  optics  books  discuss  holography.  In 
particular  the  optics  books  by  Young,  Mathieu,  and  Webb  each  provide  an  excellent 
discussion  on  the  topic  of  holography  (27, 19,  25:120-149,140-141,160-164). 

A. 2  General  Discussion 

The  term  holography  is  derived  from  the  Greek  bXo^  =  whole  or  complete  and 
'ypaipeiu  =  to  write  or  record.  TDgether  they  mean  whole  record.  Holography  can  be 
described  as  the  process  by  which  a  whole  or  complete  record  of  an  object  is  recorded. 
Usually  this  recording  is  stored  on  a  photographic  emulsion  which  is  then  referred  to 
a  hologram. 

For  a  complete  record  of  on  object  to  be  recorded  it  is  necessary  to  record  both 
the  phase  and  the  amplitude  of  the  light  which  is  scattered  by  the  object.  Since  a 
photographic  emulsion  can  record  only  intensity  it  is  necessary  to  somehow  encode  the 
phase  information.  This  encoding  is  accomplished  by  interfering  the  light  diffracted 
by  the  object  with  a  coherent  reference  beam.  When  these  two  coherent  light  sources 
interfere,  alternating  bright  and  dark  fringes  called  interference  maxima  and  minima 
are  formed. 
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Zeroth  order  maximum,  c  =  0 

First  order  maximum,  ?  =  A 
Second  order  maximum,  <;  =  2A 


Figure  11.  Young’s  double  slit  experiment 

The  classic  example  of  this  interference  is  the  double  slit  experiment  first 
performed  by  Thomas  Young  in  1802  (9:63*66).  The  experiment  consists  of  using  two 
parallel  slits  and  allowing  light  to  pass  through  the  slits  and  onto  a  screen.  The  light 
from  the  slits  will  then  produce  interference  maxima  and  minima.  Yo\mg  found  that 
the  interference  maxima  and  minima  occtirred  in  orders.  The  zeroth  order  maximum 
was  foimd  in  the  center  of  the  pattern  where  the  paths  from  the  two  slits  are  equal  in 
length  and  as  a  result,  the  path  difference,  is  zero.  The  first  order  maxima  is  found 
at  the  point  were  the  path  difference  is  A.  The  second  and  following  order  maximas 
are  foimd  where  the  path  difference  is  an  integer  multiple  of  A.  Figure  11  provides 
an  illustration  of  this.  In  general,  maxima  occur  when 

?  =  mA,  (16) 

where  m  =  0, 1,2,3, ...  and  is  referred  to  as  the  order  of  interference.  By 
convention  lA  is  considered  to  be  one  trip  around  the  unit  circle.  That  is,  the  phase 
angle  changes  from  0  to  27r  radians.  Thus,  interference  maxima  will  exist  whenever 
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the  phase  angle  difference,  6,  of  the  two  coherent  sources  is 

^  =  0,2ir,47r,6ir,...  (17) 

Conversely,  interference  minima  will  exist  whenever  the  phase  angle  difference  of  the 
two  coherent  sources  is 

S  =  Tr,3Tr,5TT, ...  (18) 

It  is  the  principle  of  interference  maxima  and  minima  described  above  that  is 
used  to  encode  the  required  phase  information  in  a  hologram.  The  fringe  pattern  of 
the  hologram  is  formed  when  two  coherent  light  sources  interfere  at  the  hologram 
plane.  One  is  the  reference  beam,  which  is  assumed  to  be  a  plane  wave.  In  other 
words,  for  all  points  (x,y)  in  the  hologram  plane  the  phase  of  the  reference  beam  is 
the  same.  There  is  no  variation.  Since  the  reference  beam  position  is  known,  and  it 
is  providing  a  plane  wave,  the  phase  is  known  a  priori.  The  second  light  source  is  the 
object  of  which  a  hologram  is  to  be  made.  Light  from  the  reference  beam  is  split  into 
two  beams.  One  is  allowed  to  proceed  directly  to  the  hologram  plane  as  discussed 
above.  The  other  is  directed  toward  the  object.  Part  of  this  light  reflects  off  of  the 
object  and  travels  to  the  hologram  plane  where  it  arrives  with  an  unknown  phase. 

It  is  at  this  point  the  phase  information  needed  to  recreate  the  image  of  the  object 
is  encoded.  Interference  maxima  will  be  formed  at  the  points  where  the  reference 
beam  and  the  object  beam  are  exactly  in  phase.  These  points  form  the  bright  fringes 
discussed  above.  Interference  minima  will  be  formed  where  the  reference  beam  and 
the  object  beam  are  exactly  ir  radians  out  of  phase.  These  points  form  the  dark 
fringes.  Thus,  since  the  reference  beam  has  a  known  phase  at  the  hologram  plane, 
the  bright  fringes  correspond  to  locations  where  the  object  beam  had  exactly  the  same 
phase.  Dark  fringes  correspond  to  locations  where  the  phase  of  the  object  beam  was 
exactly  out  of  phase  with  the  reference  beam.  Since  the  reference  beam  phase  is 
known  a  priori,  a  record  now  exist  which  indicates  the  exact  phase  of  the  light  from 
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the  object. 


A.3  Mathematical  Theory  of  Holography 

A  hologram  is  a  recording  of  an  interference  pattern.  The  interference  occurs 
between  the  amplitude  of  two  waves  Uo  and  Ur.  Uo  is  the  object  wave  and  Ur  is 
the  reference  wave.  As  Uo  and  Ur  superimpose  the  interference  fringes  are  formed 
and  recorded  on  the  photographic  emulsion.  Since  the  formation  of  these  fringes  is 
dependent  on  the  phase  and  amplitude  of  Uo  and  Ur,  it  is  crucial  that  Uo  and  Ur  be 
kept  in  a  constant  phase  relationship  during  the  recording.  That  is,  they  must  be 
coherent.  As  they  arrive  at  the  hologram  plane  Uo  and  Ur  combine  to  give  Uo  +  Ur. 
However,  the  photographic  emulsion  can  record  intensity  only.  The  intensity  at  the 
hologram  plane  is  thus  given  by 

1  =  {Uo  +  Urf 

=  iUo  +  Ur){UoAUry 

=  \Uoy+\Ur\^  +  UoU:  +  U:Ur  (19) 

where  *  indicates  the  complex  conjugate. 

When  the  photographic  emulsion  is  exposed  to  the  intensity  /,  a  transmittance 
function 

T  OC  =  (lUof  +  Wr^  +  UoU;  +  U:Ur)-^^^  (20) 

where  C  is  a  measure  of  the  emulsion  response  time. 

During  reconstruction  a  wave  identical  to  Ur  is  made  incident  on  the  hologram 
producing 

Urec  =  UrWo^  +  UrWr^  +  UrUoU;  +  UrUlUr  (21) 

where  Urtc  is  the  reconstructed  wave. 

The  first  term  of  Urec  passes  straight  through  and  is  unimportant.  The  second 


37 


term  contributes  noise  to  the  scene  and  is  not  important.  The  third  term,  UrUoU‘, 
however  produces  the  object  wave  Uq.  The  fourth  term  produces  a  second,  virtual 
image. 

A.4  Calculating  Uo 

When  producing  a  CGH  the  problem  to  be  solved  is  that  of  calculating  the 
term  Uo  in  the  equation  discussed  in  the  previous  section.  Uo  is  the  field  present 
at  the  hologram  plane  due  to  the  light  reflected  by  the  object.  In  reality  this  field 
is  continuous  across  the  hologram  plane.  Because  there  exist  no  recording  medium 
with  infinite  resolution,  and  because  the  field  is  to  be  calculated  digitally,  the  problem 
must  be  reduced  to  one  of  calculating  the  value  of  the  field  at  some  finite  number  of 
points  on  the  hologram  plane.  This  is  done  by  using  the  following  equation: 

,,,y  fexp[i2n\-'{(X.-X^f  +  (Y,-Y,r  +  (Z.-Z^m 

^  ((X. - +  (K, - +  (z, - zo*)'/''  '' 

where  X  is  the  wavelength  of  light,  Xo,  Yo,  and  Zo  are  the  coordinates  of  a  point  on 
the  object,  and  Xh,  Yh  and  Zh  are  the  coordinates  of  a  point  on  the  hologram  plane. 
H{Xh,  Yh)  is  the  value  of  the  field  at  the  point  {X,  Y)  in  the  hologram  plane.  The  term 

{{Xo  -  Xhf  +  {Yo  -  Yhf  +  {Zo  -  Zhf)^'^  (22) 

is  a  measure  of  distance  between  the  current  point  on  the  object  and  the  current  point 
on  the  hologram  plane. 

Direct  computer  implementation  of  equation  1,  in  this  form,  is  a  problem.  The 
problem  arises  because  for  this  equation  to  be  valid,  X  and  the  distance  separating 
the  object  and  the  hologram  plane  must  be  measured  in  the  same  units.  A  typical 
measure  for  X  is  .5  microns  (5  x  10"^  meters).  Thus,  if  the  object  and  the  hologram 
plane  are  at  some  point  separated  by  six  inches  in  the  X  direction,  six  inches  in  the  Y 
direction  and  six  inches  in  the  Z  direction  the  equation  to  calculate  the  partial  field 
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at  this  point  on  the  hologram  plane  would  be 


exp[i2TX-\{25, 400)2  +  (25, 400)^  +  (25, 400)2)  V2] 
((25, 400)2  +  (25, 400)2  +  (25, 400)2)  V2 
exp[i2Tr\-\l,  935, 480, 000)^/2] 

(1,935,480,000) 
eip[i27rA~^(43, 994.1)] 

43,994.1 


,i552, 846.2 


43,994.1 


No  digital  computer  available  has  the  capability  to  accurately  evaluate  e‘®®2.846.2 
Therefore,  the  form  of  equation  1  must  be  altered  to  allow  evaluation  of  the  field 
Uo  at  the  hologram  plane.  Euler’s  theorem  states  that  an  equivalent  form  of  is 
cos9  +  isin9.  In  this  case, 

9  =  2t\-\{Xo  -  X^f  +  (n  -  n)"  +  {Zo  -  (24) 


This  results  in  equation  1  becoming 


H{Xk,Yh) 


A.  5  Implementation  Note 


cos9  -h  isin9 


((x„  -  Xnf  +  (i;  -  n)2  +  {z,  -  zo2)V2 


Computer  implementation  of  equations  24  and  25  can  be  susceptible  to  numerical 
accuracy  problems.  Two  things  in  particular  must  be  avoided.  First,  one  must  insure 
that  9  does  not  become  so  large  that  the  sin  and  cos  functions  no  longer  ret;im  accnirate 
results.  Second,  simple  numeric  overflow  must  be  avoided  when  the  distances  become 
large. 
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Appendix  B.  Fresnel  Diffraction 


B.  1  Introduction 

The  equation  used  in  this  thesis  effort  is  a  general  form  for  scalar  diffraction 
theory.  It  is  not,  however,  completely  general.  Furthermore,  no  ciurent  method  for 
making  CGHs  uses  a  completely  general  form  of  the  equation.  In  all  cases,  some  form 
of  Fresnel  diffraction  is  used  to  calculate  the  intensity  values. 

This  appendix  contains  a  brief  discussion  of  the  steps  that  must  be  taken  to 
transform  the  general  form  of  the  equation  into  the  more  restrictive  case  of  Fresnel 
diffraction.  These  steps  are  taken  in  order  to  reduce  the  equation  to  a  form  that 
is  more  easily  manipulated.  The  book  Introduction  to  Fourier  Optics  by  Goodman 
provides  more  detail  on  the  subject  (ll)[57-74]. 

B.2  Approximation 

For  the  purposes  of  this  thesis  effort,  diffraction  theory  was  used  as  follows. 
Given  the  properties  of  light  at  a  given  point  P\,  determine  the  properties  of  that  light 
at  the  point  Pq,  Figure  12  illustrates  this  diffraction  geometry.  The  general  form  of 
the  equation  used  to  solve  this  problem  is 

U{xQ,yQ)  =  j  J  h{xo,yo‘,xi,yi)U{xi,yi)dxidyi  (26) 

where  f/(io,yo)  is  the  field  present  at  the  point  Po  =  (xo,yo)  as  a  result  of  the  field 
originating  at  the  point  Pi  =  (ii,  yi)  and  where 

hi^Q,yo]xi,yi)  =  -:^^^^^^^^oi)co3(n,foi)  (27) 

roi 

The  first  approximation  is  based  on  the  following  assumptions,  which  refer  to 
figure  12.  First,  we  must  assume  the  distance  z  between  the  aperture  and  the 
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yi  yo 


Figure  12.  Diffraction  Geometry  [57]goodman 

observation  plane  is  significantly  greater  than  the  size  of  the  aperture.  Furthermore, 
we  must  assume  that  from  the  observation  plane  only  a  finite  area  around  the  z  axis 
is  of  interest.  Finally,  we  must  assume  the  distance  z  is  significantly  greater  than  the 
area  of  the  region  of  interest.  For  production  of  3-D  CGHs  this  means  the  distance 
separating  the  hologram  plane  and  the  object  must  be  significantly  greater  than  the 
size  of  the  object.  If  these  assumptions  are  true  the  obliquity  factor, 

co3(n,foi) 

is  approximately  one.  This  is  the  only  approximation  made  to  the  equation  during 
tiiis  thesis  effort.  The  rest  of  the  equation  is  calculated. 

Under  the  same  conditions,  the  quantity  roi  in  the  denominator  of  equation  27 
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is  approximately  z.  This  changes  the  weighting  function  to 


H^o>yo'y^i,y\)  =  j^exp{jkroi)  (28) 

Equation  28  can  be  further  simplified  by  approximating  the  quantity  roi  in  the 
exponent.  The  distance  roi  is  given  by 


roi  = 


This  equation  may  be  approximated  by  use  of  binomial  expansion  of  the  square  root, 
producing 

yrT6  =  1  + \b\<i  (31) 

Z  o 

If  the  square  root  in  equation  30  is  assumed  to  be  adequately  approximated  by  the 
first  two  terms  of  equation  31,  the  weighting  factor  can  be  rewritten  as 

h{xo,  yo;  xi,  yi)  =  ~  +  (j'o  “  (^2) 

This  approximation  is  known  as  the  Fresnel  approximation.  When  all  of  the  above 
assiimptions  are  valid,  the  observation  plane  is  said  to  be  in  the  region  of  Fresnel 
diffraction.  Again,  the  accuracy  of  these  approximations  is  governed  by  the  size  of 
the  aperture,  the  size  of  the  viewing  region,  and  the  distance  z.  One  test  for  accuracy 
that  can  be  used  is  to  require  the  maximum  phase  change  contributed  by  the  third 
term  of  equation  31  to  be  significantly  less  than  one  radian.  This  condition  is  met 
when 

+  {yo  -  yif  (33) 


\[z^  +  (xo  -  xi)^  +  (yo  -  yi)^ 


(29) 
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If  the  Fresnel  approximation  is  valid,  equation  26  can  be  rewntten  such  that 

U{xo,yo)  =  A  j  j  U{xi,yi)exp\j^{x\  +  yl)]exp[-j2TrXz{xQXi  +  yoyi)]dxidyi  (34) 

where  A  =  ^^^^^exp[i^{xl  -f  yl)].  In  this  form  the  function  U{xq,  yo)  can  be  found  by 
taking  the  Fourier  transform  of 

U{xuyi)exp\j{^){x\  ^  yl)\ 

and  then  multiplying  the  results  of  the  Fourier  transform  by  the  amplitude  and  phase 
factors. 


Appendix  C.  Definitions 


Amplitude(A)  -  is  the  height  of  a  wave  above  the  wave’s  average  position. 
Amplitude  may  range  in  value  between  +A  and  -A.  See  figure  13. 

Intensity(I)  -  is  the  amount  of  energy  that  flows  per  unit  time  across  a  unit  area. 
The  area  is  defined  to  be  perpendicular  to  the  direction  of  the  flow.  Intensity  is  equal 
to  amplitude  squared,  I  =  A^. 

Wavelength(A)  -  is  the  distance  required  for  light  to  travel  one  full  cycle.  A 
full  cycle  can  be  measured  from  peak  to  peak  or  from  any  point  on  one  wave  to  an 
equivalent  point  on  the  next  wave.  See  figure  13. 
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-A 


Figure  13.  Wavelength  and  Amplitude 
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